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I.  INTRODUCTION 


Following  References  1  and  2,  we  propose  to  generalize  well-known  transport 
equations  for  turbulent  flows  so  that  they  are  assumed  to  apply  to  two- 
point  correlation  functions.  We  use  analogy  with  their  one-point  counter¬ 
parts  to  which  they  are  required  to  reduce  as  the  separation  of  the  two 
points  is  made  to  vanish. 


In  order  to  illustrate  the  development  of  this  program,  we  concentrate  in 
this  report  on  homogeneous  turbulence  at  high  Reynolds  numbers  so  that  the 
two  one-point  models  needed  are 


the 


dissipation  model 


.  o .  i)  (i) 


v 

and 


the  redistribution  model  t 


•  V  6i  j 


,(csl) 


(2) 


Both  models  have  been  used  successfully  for  such  a  long  time  that  it  is 
difficult  to  assign  rigidly  their  originators.  They  have  acquired  special 
importance  in  the  recent  developments  of  turbulent  transport  calculations 
(Ref.  3).  These  models  have  rather  obvious  generalizations  to  two-point  '  \ 
correlations  (they  will  be  given  in  the  next  section).  If  such  generaliza¬ 
tions  are  made,  then  the  rate  equation  for  R^j  =  <u( (x)uj(y)>  can  be 


closed  and,  in  principle,  solved,  with  the  result  that  the  structure  as 
well  as  the  scale  of  the  turbulent  field  can  be  determined. 


A. 


For  homogeneous  flows 
of  Rij 


,  using  the  leading  term  in  the  moment  expansion 


to  obtain  a  closed  scale  equation  is  permissible  because  the  correlation 
length  is  infinitely  small  compared  with  the  size  of  the  system. 


1.  Sandri,  G.,  "A  New  Approach  to  the  Development  of  Scale  Equations  for 
Turbulent  Flows,"  A.R.A.P.  Report  No.  302,  April  1977. 

2.  Sandri,  G.,  "Recent  Results  Obtained  in  the  Modeling  of  Turbulent  Flows 
by  Second-Order  Closure,"  AF0SR-TR-78-0680,  February  1978. 

3.  Donaldson,  C.duP.,  "Construction  of  a  Dynamic  Model  of  the  Production 
of  Atmospheric  Turbulence  and  the  Dispersal  of  Atmospheric  Pollutants," 
DloikAhop  on  HicAomeXzoKology  (D. A. Haugen,  ed.),  American  Meteorological 
Society,  Boston  (1973),  pp.  313-392. 


B.  Experimental  information  in  available  on  several  variables  for  both 
isotropic  grid  turbulence  and  homogeneous  shear  turbulence  (Refs.  4 
and  5). 

With  only  two  adjustable  constants,  the  model  covers  qualitatively  both 
types  of  turbulent  flows.  We  shall  show  this  by  exhibiting  explicit 
analytic  solutions  with  several  of  the  desired  features.  In  particular, 
the  analytic  solutions  for  the  homogeneous  turbulence  models  show  the 
presence  of  two  distinct  time  scales  which  characterize,  respectively,  the 
rapid  settling  of  the  tensor  character  of  the  flow  to  an  asymptotic  state 
and  the  slower  development  of  the  energy  and  mean  scale.  It  is  found  that 
for  both  grid  and  homogeneous  shear  turbulence,  the  ratio  of  the  two  scales 
is  about  ten. 

The  fast  time  is  the  redistribution  time,  A/q  ,  while  the  slower  one  is 
bA/q  (dissipation  scale)  for  grid  turbulence  and  v'A/q  (merging  scale)  for 
shear  flows  (v1  is  defined  in  the  next  section)  after  an  initial  transient. 
This  feature  of  the  model  solutions  seems  to  be  well  reflected  in  the  data. 

In  this  report,  we  will  construct  the  general  structure  equations  (Appendix 
A).  We  shall  also  obtain  a  first-order  test  of  our  model  by  exploiting 
the  following  result  for  homogeneous  turbulence: 

The  equation  for  R-j j  is  wholly  determined  from  three  requirements: 

(i)  that  it  should  yield  the  observed  transport  for  u-jUj  ; 

(ii)  that  it  should  satisfy  the  continuity  equation; 

(iii)  that  it  should  yield  the  correct  limit  for  isotropic  turbulence. 


We  now  consider  briefly  an  analogy  between  the  models  (1)  and  (2)  and  the 
Newton-Fourier  heat  equation.  We  may  think  of  the  Newton- Fourier  equation, 

=  -K  %  <3> 

which  gives  the  heat  flow  vector  q-j  in  terms  of  the  temperature  gradient, 
as  a  phenomenological  law  (or,  more  incisively,  "model")  which  allows  us  to 
close  the  heat  equation  and  hence  gives  us  a  chance  to  solve  it.  As  a 
model,  (3)  is  subject  to  restrictions.  It  is,  however,  tensorial  and  there¬ 
fore  independent  of  geometry,  hence  valid  in  any  coordinate  system.  Fourier 
luckily  relied  (intuitively)  on  the  tensor  nature  of  (3)  and  overlooked  the 
restrictions:  (i)  K  is  temperature-dependent  even  for  the  simplest 
material  (inert  gases);  (ii)  K  is  a  tensor  for  any  nonsimple  material; 
(iii)  in  the  presence  of  electricity,  qi  requires  altogether  a  new  term 
(the  thermoelectric  effect  of  Thomson). 

4.  Harris,  V.G.,  Graham,  J.A.H.  and  Corrsin,  S.,  "Further  Experiments  in 
Nearly  Homogeneous  Turbulent  Shear  Flow,"  J.  Fluid  Mech.  81,  1977,  pp. 
657-687.  Corrigendum,  J.  Fluid  Mech.  86,  1978,  pp.  795-755”. 

5.  Corrsin,  S.  and  ((oilman,  W.,  "Preliminary  Report  on  Suddenly  Sheared 
Cellular  Motion  as  a  Qualitative  Model  of  Homogeneous  Turbulent  Shear 
Flow,"  Proc.  SQUID  Symp.  on  Turbulence  in  Internal  Flows,  pp.  11-33 
(S.N.B.Murthy,  ed.).  Hemisphere  Publishing  Co.,  1977. 
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Of  course,  analogous  remarks  may  well  apply  to  (1)  and  (2).  It  seems  to 
us  that  Fourier's  work  encourages  the  view  that  a  good  treatment  of  the 
simplest  model  is  desirable. 


■  v"  ■* 


II.  DERIVATION  OF  THE  CLOSED  EQUATIONS  FOR  THE  REYNOLDS 
AND  SCALE  TENSORS 


We  start  with  the  Navier-Stokes  equations  for  the  velocity  field  and 
the  kinematic  pressure  p  (pressure/density) 


at 


k  ax. 


%  ■ 

(4) 

au. 

*t*° 

(5) 

The  velocity  and  pressure  are  decomposed, foil  owing  Reynolds,  into  mean  and 
fluctuating  parts 


u  •  +  u  • 
l  i 

(6) 

p  +  p* 

(7) 

Substituting  (6)  and  (7)  into  (4)  and  (5),  one  obtains  separate  equations 
for  the  mean  and  fluctuations.  After  standard  operations,  these  equations 
can  be  cast  into  the  form  of  equations  for  the  correlation  tensor  R-jj 

Rij(x,y)  =  <u!(x)ul(y)>  (8) 


namely. 


at 


(9) 


4 


mmmm 


and 


3^  =0 

(10) 

3Rii 

*f=° 

(11) 

It  is  convenient  to  recast  these  equations  by  introducing  "centroid"  and 
"relative"  variables  as  follows: 


x  .2-+* 
c  2 


r  =  y  -  x 


(12) 

(13) 


Using  the  chain  rule  for  differentiation,  the  equations  for  R-jj  become 


3t 


Uk(x)  +  u.  (y)  a  ,  \  s 

2  3x^  +  (uk(y)  "  uk(x))ar^ 


Rij  + 


3u.  3u. 

1  3  [<»J(«K(x)«i(y)>  +  <“i(x)u^(y)uj(y)>] 


ck 


*k 


^u! (x)u£(x)uj(y)^>  -  <yu!(x)u£(y)uj(y))>] 


lj-7  <^P'(x)uj(y)^>  +  f^—  ^u-(x)p'(y)^> 


ci 


cj 

/  3uHy)\  /3u!  (x)  \ 

\p,(x)^r) +  \^rp,(y)/J 

•  (y)\ 

ck  / 


+  VVC  Rij  *  2V\3X 


3u! (x)  9ul(y)' 


(14) 
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1 


*■"-  -vro  -- 


3R^  3Ri . 
3*ci  '  2  ^ 


=  _2  d-L 

3xcj  3rj 


where  the  arguments  in  the  correlations  are  given  by 


x  =  xc  "  2 


For  homogeneous  turbulence,  the  derivatives  with  respect  to  the  centroid 
vanish  when  acting  on  any  correlation  (but  not,  in  general,  on  the  mean 
velocity)  or  on  a  fluctuation.  We  then  have,  for  homogeneous  turbulence. 


3Rii  T.  .  ] 

=  -  M*)  -  M*)  ^r1  -  R 


3u^  3u- 

ik^  +  ^"kj  + 


fCui(x)uk(x)Uj(y)^  -  <u:(x)u'(y)uj(y)>] 


+  (  P‘(x) 


3u!(y) 


’au'.(x)  3uj(y) 
3xck  "  KT 


3u!(x) 


p'(y) 


The  terms  in  Eq.  (19)  that  prevent  closure  of  the  equation  are  the  last 
three.  In  order  to  achieve  closure,  we  introduce  generalized  transport 
models  as  follows: 


Dissipation  equation 


-2v 


Sxck  Sxck  /  HaJs  'W,'> 


(22) 


Intercomponent  rearrangement 


P 1  (x) 


3uHy)  3u!(x) 


9xci 


3x 


p’(y) 


cj 


1 

A 


(23) 


where  "L  •  is  assumed  to  satisfy 

*  J 


Ti j(r  -  0)  -  0 


Ti j(r) 


d3r  ■  0 


(24) 

(25) 


These  properties  are  utilized  in  Section  III  to  derive  the  Reynolds  stress 
and  scale  tensor  equations.  For  a  full  determination  of  eddy  structure, 
additional  properties  are  needed  which  will  be  discussed  in  Appendix  A.  In 
addition  to  the  two  generalizations  given  above,  we  need  a  model  for  the 
spatially  homogeneous  part  of  the  triple  velocity  correlation  which 
represents  the  nonlinear  effects  of  local  turbulent  convection.  These  non¬ 
linear  effects  correspond  to  either  eddy  break-up  (cascading,  when  wave 
vectors  add)  or  merging  of  eddies  (when  the  wave  vectors  subtract).  We 
thus  call  the  model 

Eddy  size  rearrangement 


(y)^>  -  <(u!(x)u£(y)uj(y))> 


where  N. *  satisfies 

*  J 


N..(r  «  0)  =  ujuj 


(27) 

(28) 
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where  the  Dirac  function  of  the  magnitude  of  r  is  related  to  6(r)  by 

6(r)  *  A1M1  (33) 

4tfr 

The  approximation  (33),  as  we  remarked  in  the  Introduction,  is  best  justi¬ 
fied  for  homogeneous  turbulence  because  in  this  case  the  spatial  scale  of 
the  mean  flow  is  infinite. 

We  now  apply  to  (29)  the  operation  f d^r/4nr^  and,  using  (25)  and  (28),  we 
find  J 


-q‘ 


3u, 


Aik  3x 


ck 


a 

A 


Aj  -  5  6lA 


kk 


21Aj 


(34) 


Convection  of  A-jj  does  not  occur  for  homogeneous  turbulence.  To  see  that 
this  fact  is  a  consequence  of  the  moment  approximation  (32),  we  note  that 
(33)  gives 


where  we  performed  an  integration  by  parts  and  used  continuity.  We  obtain 
the  final  form  of  the  equation  for  the  scale  tensor  by  substituting  into 
(34) 


3_ 

at 


a„2 

A.  .  29- 

ij  at 


(35) 


9 


(36) 


tX 


-  9uk  a3 

-2ui“k9F7-  2b^r 

Cl 


This  latter,  the  energy  equation,  is  obtained  by  contracting  (30).  The 
result  is 


^ A. 

at 


au. 


3ui 


ik  3x 


ck 


ax 


ck 


kj 


+  Aij 


at  +  *(2b  +  v'> 


.  £ 

A 


2bq5ij 


(37) 


Note  that  Eq.  (37)  is  obtained  by  dividing  by  q^  and,  therefore,  it  should 
not  be  used  in  the  absence  of  turbulence. 

To  close  the  pair  (30)  and  (37),  we  choose 

*-Kk  (38> 


The  reason  for  this  choice  was  discussed  in  a  previous  report  (Ref.  1). 
Contraction  of  (37)  and  use  of  (38)  yields  the  equation  for  the  mean  scale 


3A  _ 

at 


2  3u 

3  Ai k  ax 


i  2ukui  9llk 
1  +  A  +  v'q 


ck 


.2  ax 


cl 


(39) 


It  is  interesting  to  note  that,  in  this  model,  the  coefficient  of  the 
"production"  term  is  not  a  universal  constant  and  receives  generally  compet¬ 
ing  contributions  from  the  Reynolds  stress  and  from  the  tensor  scale. 
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III.  SOLUTIONS  OF  THE  COUPLED  EQUATIONS 
FOR  STRESS  AND  SCALE  TENSOR 


In  this  section  we  give  two  analytic  solutions  to  the  coupled  equations  for 
the  stress  tensor  and  scale  tensor  equations.  We  use  subsections  to 
separate  the  different  calculations. 

A.  Equations  in  Standard  Coordinates 

The  centroid  vector  and  mean  velocity  vector  are  taken  to  have  components 


(x,y,z)  ,  ( V (y )  ,0,0) 

with  3U/9y  =  U‘  =  constant.  The  relevant  components  of  the  Reynolds  stress 
equations  are  obtained  from  Eqs.  (30)  and  (36).  We  drop  primes  on  the  __ 
fluctuations  and  give  a  form  useful  for  numerical  integration  in  which  u7 
and  are  calculated  from  1 


~2 


~2  72 


U1  =  ^  -  u2  ‘  u3 


(40) 


A11  3A  A22  "  A33 


(41) 


The  other  relevant  components  of  the  stress  and  energy  equations  are 


8t  3 


9t  ulu2 


la 

dt 


i (1  ■ 

2b)  ar'  a  “2 

(42) 

3“- 

(43) 

1 

C 

c: 

’  '  A  ul“2 

(44) 

ulu2 

q 

2 

U'  -  b  3- 
A 

(45) 

For  the  tensor  scale  components,  we  obtain,  using  Eqs.  (37)  and  (39) 


22 


j  A22  +  (1  -  2b)q 


(46) 


11 


(47) 


i_  a 

at  a33 


TA33  +  •  2b>« 


at  A12  =  "  T  A12  '  U'A22 


(48) 


where 


3  uiup  2 

ft  A  '  2  ~^U'A  -  f  A12U'  ♦•'I 


1  u. u9  . 
y  =  -2  U'  +  J  (1  -  2b  -  v') 

q  n 


(49) 


(50) 


B.  Solution  of  the  Shearless  Equations 

Setting  U'  =0  ,  we  see  that  equations  for  q  and  A  decouple  from  the 
tensor  components.  Introducing  the  deviators 


_  l  2 

d. .  =  u.u.  -  6.  -q 

U  i  3  3 


(51) 


°ij  Aij  "  3  6ijAkk 


(52) 


and  the  time 


A 

T  =  q 


(53) 


We  have  the  set 


a  b  a  ,  v'  , 

at  q  °  "  t  q  •  at  A  “  Ta 


a j  .  b_  ,  a_n  .  l  -  2b  -  v1 

at  ij  "  "  t  dij  *  at  ij  t 


fyt=b+V’ 


IJ 


(54) 

(55) 

(56) 


Integrating  (56)  as 


t  ■  (b  +  v 1 )  ( t  -  t  )  +  =- 

^0 


(57) 
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we  see  that  q  ,  A  ,  d-jj  and  D-jj  are  suitable  powers  of  (q0/T0)t  ;  for 
example,  J 


q  =  qr 


d.  .  =  d. . 

U  iJ, 


A  =  A. 


(b  +  v1)  /  (t  -  tj  +  1 
0 


1-b/(b+v') 
-l/(b+v* ) 


(b  +  V*)  -O  (t  -  t)  +  1 
0 


-1  -v '  /  (b+v  * ) 


(b  +  v')  yT  (t  -  t0)  +  1 
0 


(58) 

(59) 

(60) 


From  grid  turbulence  data  on  q  and  A  ,  we  may  choose 


• 


We  then  see  that  for  large  times 


q  ~  Qr 


v*  =  0.075 


(61) 


(62) 

(63) 


which  shows  that  the  deviator  decays  with  a  power  about  four  times  larger 
than  the  energy.  This  substantial  difference  may  eventually  be  checked  in 
our  anisotropic  grid  -now. 

From  the  solutions  given  above,  we  can  verify  that  statistics  are  preserved 
by  the  model  equations  if  the  model  parameters  satisfy  certain  bounds.  We 
first  show  that  the  two  tensors  u^Uj  and  A-jj  are  positive  definite  from 
their  definitions.  Consider  an  arbitrary  (constant)  A^  then. 


AiuiUjAJ  =  (u-A)2  >  0  (64) 


the  equality  sign  holding  for  A  =  0  only.  Thus  u^Uj  is  a  positive 
definite  tensor.  J 

From  the  definition  of  the  scale  tensor  given  by  Eq.  (31),  using  Fourier 
transform  on  R. .  , 

*  J 


13 


tA.,  =  l^r -JS 


ari  d3k 

8iTk 


(65) 


where  the  power  spectrum  tensor  4>ij  is  positive  definite  by  Khiutchine  s 
theorem.  Thus, 


V  AiAijAj  =  JsiirAi*ijAj  - 


(66) 


Thus,  A-: -i  is  positive  definite  because  q2  is  positive  as  a  consequence 
of  (64).  J 


Using  the  solution  (61)  and  an  analogous  solution  for  A-jj  ,  we  find 


-l/(b+v' ) 


Vj(t)  =  Vj(0)( 


+  lfi  q2 
3  iJMo 


■  v2b/(b+v,)  (%  vI/(6tv,) 


A 


A. 


Vt)  =  V0)(^T 


o  /  \  o 

-(l-2b-v' )/(b+v' ) 


(67) 


*  msiS 


v ' / (b+v ’ ) 


-(l-Zb-v’JAb+v’)' 


o  / 


(68) 


We  now  multiply  (67)  by  A.Aj  when  Ai  is  an  arbitrary  vector  and  find 

\  -l/(b+v' ) 


AiVj(t)Aj  =  (A-ur(0)(  fi 


A2  _  2 


+  T% 


’  \  -2b/(b+v‘)  ,  \  -1/ (b+v ' ) 

A-  •<•) 


(69) 


From  (57)  we  see  that 


^  T  >  1 

Ao  “ 


(L  +  v'  >  0) 


(70) 
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Sufficient  for  the  left-hand  side  of  (69)  to  be  positive  is 

\  -2b/ (b+v ' )  ,  v - 1/ (b+v * ) 

)  *(<’) 

which  requires,  using  (70), 

2b  <  1  (72) 

A  similar  analysis  applies  to  Aij  ;  however,  no  further  restrictions  on 
the  parameters  is  found. 

C.  An  Exact  Solution  of  the  Equations  with  Shear 
To  obtain  a  solution  of  the  equations  with  shear,  we  let 


9  = 

V  eav  t 

A  = 

L  t**'1 

u,  ,2a,,t 

All- 

Lj  eavt 

J 

II 

U2  e2avt 

A22  = 

4  ea,,t 

-J 

II 

K,  e2avt 

A33  = 

4  eavt 

iU2  = 

u4  e2av'1 

A12  = 

L,  ea»'* 

(73) 


Substituting  these  forms  into  the  differential  equations  of  III. A,  we  find 
that  the  exponentials  cancel  and  that  an  algebraic  set  of  equations  for  the 
amplitudes  are  obtained.  It  is  possible,  with  some  algebra,  to  solve  the 
amplitude  equations  explicitly  in  terms  of  the  parameters  b  and  v'  . 

The  energy  components  are 


U1  _  1  +  6v'  +  4b 

t  'am  +  jv'T 


(74) 


2  2 

u2  .  u3  „  1  -  2b  I 

7  *7  —y— r+r^ 

9  9 


(75) 


The  scale  components  are 


A11  _  1  +  6v'  +  4b 
"A - PTTT1 - 


(76) 
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The  off-diagonal  components  are 


The  Corrsin  parameter  is 


The  ratio  of  the  two  times  is 


and  the  growth  rate,  a  ,  is 


We  make  the  following  remarks: 


1)  Several  numerical  integrations  suggest  that  any  solution  that  initially 
statistically  realizable  will  remain  so  and  will  asymptote  to  the 
exact  solutions  given. 

2)  If  the  scalar  scale  equation 

|=c^u^U'  +v‘q  (c  »  0.35)  (86) 

is  adopted  instead  of  the  tensor  equation,  an  exponential  solution 
exists  and  has  the  same  qualitative  features  insofar  as  scale  and 
energy  are  concerned  (of  course,  no  scale  directivity  results).  This 
indicates  that  the  models  for  scale,  which  were  not  designed  to  fit 
homogeneous  shear  data,  are  quite  stable. 

3)  The  equations  for  Rjj  itself  (or  the  spectrum  <j>ij  ,  see  Appendix  A) 
have  the  same  type  of J convective  equilibrium  solutions. 


IV.  CONCLUSIONS 


Our  analysis  of  a  model  of  homogeneous  turbulence,  deduced  from  simple 
assumptions  made  on  the  two-point  tensor  R-jj  ,  has  brought  out  four  main 
conclusions. 

1.  The  model  contains  only  two  adjustable  parameters  which  we  chose  to 
fix  from  grid  turbulence  data.  Good  qualitative  agreement  with  homo¬ 
geneous  shear  flow  results. 

2.  The  model  implies  the  existence  of  two  distinct  time  scales  whirh  are 
separated  by  a  factor  of  about  10.  They  appear  after  an  initial 
transient  phase  has  died  out  during  which  q/A  becomes  approximately 
equal  to  U‘  .  On  the  fast  scale,  A/q  ~  (U * )~ 1 ,  the  normalized 
deviator 


locks  into  a  constant  value  indicated  by  a  Corrsin  parameter 


0.5 


or  by  a  Bradshaw  number 


On  the  slower  scale,  v'(A/q)  -  0.07(A/q)  ,  the  energy  components  and 
the  mean  scale  grow  exponentially.  This  solution  can  be  thought  of  as 
a  superequilibrium  with  convection.  A  qualitative  physical  picture  is 
as  follows.  During  the  initial  transient,  when  (qo/Ao)«U'  ,  the 
shear  brings  the  sudden  distortion-like  turbulence  up  to  convective 
equilibrium  (q/A  *  U * ) ,  while  when  (qo/Ao)  »  U‘  ,  the  turbulence  decays 
grid-like  to  convective  equilibrium.  Then  a  merging  mechanism  takes 
over  (a  multi-layer  Brown-Roshko  effect)  so  that  eddies  fold  with  each 
other,  making  larger  ones  indefinitely  (as  long  as  the  imposed  shear 
provides  the  energy  to  sustain  the  process).  Once  the  merging  process 
takes  over,  the  eddy  structure  remains  fixed  and  exhibits  highly  direc¬ 
tional  integral  scales  (see  4  below). 

3.  In  the  model ,  the  two  transverse  energy  components  do  not  separate 
while  the  data  Indicate  that  such  a  separation  occurs.  The  splitting 
can,  however,  be  brought  into  the  model  by  assuming  a  tensor  coefficient 
In  the  redistribution  equation.  We  have  taken  this  point  to  constitute 
a  refinement  at  this  stage  of  analysis. 
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APPENDIX  A.  SPECTRAL  EQUATION 


The  equation  for  the  spectral  tensor  is 

*  J 


» 


*  •  U'(*i25jl  ♦  S2i*2j>  -  1  [V*)  ♦  V-*>] 

*  t  [k,PjW  -  Kj?1(.J)]  -  2vk2*fj 


where 


VK)  =  Ka  Siaj(r)  e 


■iK*r  dr 


PjU)  -  J*<P'(x)uj(y)>e-iK* 


r  dr 


Continuity  requires  the  Poisson  equation 


2  P,(k)  -  Pj(-K) 
iic  -J - g— J - 


( A  - 1 ) 


‘  2U  Va  * V‘>  +  V'*1]  (A-2) 


The  following  simple  model  gives  a  tolerable  picture  for  the  energy 
spectrum  of  grid  turbulence: 


‘  a  j  * 6  k  fe  +u  *  2*t j  * c  (fc2  i  * 4k  If  »(j  *  z*i j) 

(A. 3) 


A  s  0.075  ,  B  =  0.615  ,  C  =  0.27 


which  have  been  adjusted  to  give  exactly  a  5/3  law  for  equilibrium. 

Even  though  a  complete  structure  equation  has  not  been  determined  for  the 
sheared  turbulence,  it  is  seen  that  any  closure  scheme  such  as  that  given 
by  Eq.  {A. 3)  will  have  the  exponential  solution  of  the  form  given  in 
Section  III. 


